Astronomical Techniques and Instruments, 2023 
Article Open Access 


Determining the observation epochs of star catalogs from 
ancient China using the Generalized Hough Transform method 


Boliang He!2 Ô, Yongheng Zhao!2* © 


‘National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China 


2University of Chinese Academy of Sciences, Beijing 101408, China 


* 
Correspondence: yzhao@nao.cas.cn 


Received: January 3, 2024; Accepted: January 20, 2024; Published Online: March 8, 2024; https://doi.org/10. 
61977/ati2024013 


© 2023 Editorial Office of Astronomical Techniques and Instruments, Yunnan Observatories, Chinese Academy of 
Sciences. This is an open access article under the CC BY 4.0 license (http://creativecommons.org/licenses/by/4.0/) 


Citation: He, B. L. and ZHAO, Y. H.. 2023. Determining the observation epochs of star catalogs from ancient China 
using the Generalized Hough Transform method. Astronomical Techniques and Instruments, littps://doi.org/10. 
61977/ati2024013. 


Abstract: Ancient China recorded a wealth of astronomical observations, notably distinguished by the inclusion of 
empirical measurements of stellar observations. However, determining the precise observational epochs for these 
datasets poses a formidable challenge. This study employs the Generalized Hough Transform methodology to analyze 
two distinct sets of observational data originating from the Song and Yuan dynasties, allowing accurate estimation of 
the epochs of these stellar observations. This research introduces a novel and systematic approach, offering a scholarly 
perspective for the analysis of additional datasets within the domain of ancient astronomical catalogs in future 
investigations. 
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1. INTRODUCTION 


Records from ancient China include a wealth of astro- 
nomical observations, prominently featuring a valuable com- 
pilation of meticulously measured stellar observations. Con- 
temporary scholars have analyzed these data extensively, 
but the observational epochs remain elusive. Conven- 
tional methodologies use the precession correction 
approach, regressing observational data with precession 
datal!l. Fourier analysis methods have also been applied 
to estimate epochs by examining stellar data to scrutinize 
declination errors!?]. 

Previous methodologies for epoch determination 
required highly accurate correspondences between ancient 
and contemporary stars, with their precise positional coordi- 
nates. In this study, we propose the implementation of a 
novel algorithm using the Generalized Hough Transform 


to enhance the precision of calculations], Unlike its prede- 
cessors, this technique operates through probabilistic calcu- 
lations, alleviating the requirement for absolute accuracy 
in the data, with substantial correctness in the majority of 
the data deemed sufficient. In the computational process, 
only the departure from the pole degree data of the stars 
is essential, differentiating it from other methods that man- 
date not only the departure from qu ji du & $k FE (dis- 
tance to the pole) but also ru xiu du A fñ FE (distance to 
the west boundary of xiu). 


Fig. 1. Illustration of an ancient Chinese equatorial coordinate 
system. 


This study employs the Generalized Hough Trans- 
form to analyze ancient stellar observations, aiming to pre- 
cisely estimate their associated epochs. 


2. TECHNICAL METHODS 


The Hough Transform is an ingenious method that 


transforms a global curve detection problem into an effi- 
cient peak detection problem in parameter space. It accom- 
plishes this by identifying imperfect instances of objects 
with a specific type of shape through a voting process. 
This voting process occurs in a parameter space, wherein 
candidate objects are regarded as local maxima within the 
accumulator space explicitly constructed by the algorithm 
responsible for computing the Hough Transform. 


2.1. Algorithm Principle 


The records of celestial body positions in ancient Chi- 
nese astronomy are commonly expressed in terms of two 
parameters. The ru xiu du A {ii FE is the distance to the 
west boundary of xiu, astronomically equivalent to the 
right ascension of a star in an equatorial coordinate sys- 
tem. The only distinction between this and the modern sys- 
tem lies in its reference point, which is not the spring 
equinox but the determinative star (ju xing) of the xiu. 
This coordinate is always expressed as the number of 
degrees within a specific xiu. The other parameter, the qu 
jidu EWE or distance to the pole, is equivalent to decli- 
nation. Using Fig. | as an illustrative example, point X 
denotes the position of a star within a specific mansion, 
and P designates the north celestial pole. Consequently, 
PX signifies the reference meridian for this mansion. 

For a star S in this mansion, the difference in longi- 
tude 0 represents the entry degree, while the angle œ 
from the north celestial pole to point S is the departure 
degree. Ru xiu du (RXD) is denoted by 0, and qu ji du 
(QJD) is represented by œ. The relationship between the 
QJD œ and the declination (6) in the modern equatorial 
coordinate system is given by 


6=90°-w 


In contrast, in modern observational astronomy, stel- 
lar positions in the celestial sphere are usually articulated 
in terms of right ascension, a, and declination, ð, as 


S = (Qs, ôs) 


In the Q/D data, œ signifies the angular distance of a 
star from the north celestial pole. The position of the 
north celestial pole at the observational epoch in the mod- 
ern equatorial coordinate system is denoted as p = (a,6). 
Therefore, for any given star, the dot product of the two 
vectors gives 


s: p = |s||plcos w = cos w (1) 


In the (a,ô) parameter space, the north celestial pole 
for each set of data corresponds to a circle, as shown in 
Fig. 2. Nine sets of stellar observation data near the north 
celestial pole have been selected, and it is evident that the 
majority of the circles in the parameter space converge to 
a singular point. Only one circle deviates slightly from 
the convergence point. With an increased dataset, the opti- 
mal convergence point for data can be ascertained proba- 
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bilistically, serving as the candidate location for the north 
pole. 
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Fig. 2. The parameter space of stellar positions after Hough 
Transform, showing spatial location. The blue lines depict the 
representation of nine stars in the parameter space. 


The Hough Transform can be used to process and cal- 
culate all of the stellar observation data, allowing the deter- 
mination of the coordinates of the north celestial pole at 
the observational epoch in the modern reference frame. Sub- 
sequently, by incorporating the position of the north celes- 
tial pole along with precession data (refer to Fig. 3), the 
true epochs of the stellar observations can be ascertained. 
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Fig. 3. The position of the north celestial pole from 500 BCE to 
2100 CE. 


The rotational axis of Earth undergoes precession in 
an approximately 26 000-year cycle, and the varying posi- 
tions of the north celestial pole at different epochs are illus- 
trated in Fig. 4. Historical positions of the north celestial 
pole are calculated using the precession algorithm devel- 
oped by Jan Vondrákl4 51, 


Fig. 4. A schematic diagram solving the above system of 
equations, and also depicting the solutions in the parameter 
space after the Hough Transform. 


2.2. Calculation Methods and Steps 


We assume that the observational records within each 
dataset are generated over a relatively short time period, 
implying that their QJD values are grounded on the same 
zero point. These records are denoted by their coordi- 
nates in the J2000 epoch as (a; ‚ð; (i = 1, 2,..., m), and 
can be converted into Cartesian co-ordinates as (x; y; , Z; ) 
(i = 1, 2,..., m). The associated QJD data are represented 
as œw; (i= 1, 2,..., m). 

Assuming the position of the north celestial pole dur- 
ing this time period is denoted as (a,0), it can be con- 
verted to Cartesian coordinates as p (x,y, Z). 

For any two stars, A and B, each characterized by spa- 
tial positions s; (x; „Yi , Z) and s; (x; , yj Z; ), where i # j, 
the relationship is expressed as 


Si’ P =COSW; 
Sj’ Pp=COSWj (2) 
Ipl=1 


which can be represented in the Cartesian coordinate sys- 
tem as 


XiX +yiy + ZZ = COS Wi 
XjX+YjY+ZjZ=COSWj (3) 


xy er =1 


To resolve the system of equations described above, 
as illustrated in Fig. 3 in which the points on the blue 
line represent the positions of the north pole for different 
years, three potential scenarios for solutions exist: a 
unique solution, double solution, and no solution. Employ- 
ing iterative methods facilitates the calculation of intersec- 
tion solutions between any two stars. 

Considering potential errors in observational data, it 
becomes imperative to introduce a certain degree of error 
sampling during the process of solving the system of equa- 
tions. This practice proves effective in augmenting the preci- 
sion of the results. 

Owing to the precision issues in historical data observa- 
tions, it is possible to introduce a certain amount of ran- 


dom error into the calculations for modeling purposes. By 
incorporating a large amount of random data within a cer- 
tain margin of error, it is possible to enhance the results 
of data calculations. The result, observed in the phase 
space of right ascension and declination, corresponds to 
the maximum probability value near the final position of 
the north pole. 


2.3. Preprocessing of the Raw Data 


Prior to data processing and analysis, the declination 
data need to be converted. The degrees in the observa- 
tional data are represented in the traditional ancient 
degree system (gudu th J, Chinese Degree), wherein a 
full circle corresponds to 365.25° in the celestial sphere. 
Consequently, as part of the data processing, a conver- 
sion is required to align the data with the modern system, 
where a full circle corresponds to 360°, using the relation 


3. REFERENCE DATA 


During the data processing, we used the latest celes- 
tial measurement framework and modern stellar position 
data for our calculations. The calculation library 
employed in this article is the IAU’s Standards of Funda- 
mental Astronomy (SOFA) 7], which guarantees accu- 
rate results. 

The modern stellar catalog data used in our study are 
sourced from the Hipparcos Catalog!*!, the updated Hippar- 
cos Catalog from 200701, and the extended Hipparcos Cata- 
log (XHIP) from 2012[!°l, XHIP is derived from the Hippar- 
cos Catalog but integrates the latest stellar proper motion, 
parallax, and radial velocity data. This comprehensive 
dataset allows for the accurate calculation of stellar posi- 
tions across various historical epochs. The epoch for these 
catalogs is J2000.0!!!1. 

Once an approximate date for an ancient star catalog 
is established, further calculations incorporating proper 
motion and other relevant data are used to precisely deter- 
mine the positions of stars within that catalog during the 
specific historical era. This additional step serves to vali- 
date and confirm the epoch of the data under considera- 
tion. 


4. DATA PROCESSING AND ANALYSIS 


We conduct an analysis on the two sets of star cata- 
logs using the algorithm described above. 


4.1. Song Dynasty Star Catalog 


During the Northern Song Dynasty (960-1127 CE), 
there was a flourishing of society, the economy, and scien- 
tific and cultural pursuits. In the realm of astronomy and 
calendrics, astronomers actively engaged in the manufac- 
ture and calibration of astronomical instruments, undertak- 
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ing extensive observations of celestial phenomena. Their 
remarkable achievements and records stand as a testa- 
ment to the zenith of ancient Chinese technology and 
knowledge during that era. 

Among the extant historical documents from the 
Song Dynasty, notable works such as Lingtai Miyuan #i 
== Bh 91, Wenxian Tongkao Chk i 7%, Guankui Jiyao & 
Mi, Oianxiang Tongjian Y Ri, and Tianyuan Lili 
K Jt i #4 contain a wealth of stellar observation records. 
Pan Nai and othersl!? !3] meticulously organized and ana- 
lyzed these documents, including a Ming Dynasty copy of 
Lingtai Miyuan Œ 2% M4 3E housed in the National Library 
of China. Their efforts resulted in the compilation of a stel- 
lar catalog containing data for 360 stars, known as the 


Table 1. Results of the Song Dynasty star catalog processing 


Song Huangyou Star Catalog R =% Æ X. In this section, 
this catalog serves as the foundation, and the Hough trans- 
form algorithm is applied for the purpose of reprocessing 
it to ascertain its epoch. 

The calculation results are presented in Table 1, 
where P4i denotes the distance from the candidate point 
to the nearest north pole. Fig. 5 illustrates that selecting dif- 
ferent values for o yields distinct yet approximately simi- 
lar candidate outcomes. Notably, in the cases of o = 0.25 
gudu and o = 0.3 gudu, the convergence of candidate 
points appears relatively favorable. Fig. 6 specifically repre- 
sents the scenario with o = 0.25 gudu. Consequently, it is 
reasonable to infer that the epoch falls approximately 
between 1011 CE and 1016 CE. 


Data Source Verified value Lingtai Miyuan1 Lingtai Miyuan2 Wenxian Tongkao Xianglin GuankuiJiyao Tianyuan Lili 
Stars 352 342 347 263 297 340 300 
QJD error Epoch Pict Epoch Pgs Epoch Pgs Epoch Pgs Epoch Pg Epoch Pgs Epoch Pst 
(gudu) (deg) (deg) (deg) (deg) (deg) (deg) (deg) 
0.10 995 0.03 989 0.11 987 0.07 1003 0.07 998 0.07 983 0.11 1002 0.08 
0.20 1005 0.04 996 0.01 996 0.06 1015 0.06 1006 0.05 997 0.15 1012 0.07 
0.25 1011 0.04 1002 0.00 1002 0.06 1022 0.05 1012 0.04 1009 0.16 1019 0.06 
0.30 1016 0.05 1007 0.01 1008 0.07 1028 0.04 1018 0.03 1019 0.17 1026 0.05 
0.35 1019 0.05 1010 0.02 1012 0.07 1033 0.04 1025 0.03 1026 0.18 1033 0.04 
0.40 1021 0.05 1013 0.03 1016 0.07 1034 0.03 1028 0.03 1029 0.17 1036 0.04 
0.45 1022 0.06 1013 0.03 1019 0.08 1033 0.02 1032 0.03 1030 0.16 1036 0.03 
0.50 1022 0.06 1014 0.04 1020 0.08 1035 0.02 1031 0.03 1030 0.16 1036 0.03 


In the year 1011 CE, during the fourth year of 
Emperor Zhenzong's Dazhong Xiangfu XK F # ff reign in 
the Song Dynasty, as documented in Song Huiyao Jigao 
RA BEA Fig], historical records note, "In the third year 
of Dazhong Xiangfu, on the leap second of the second 
month, the Bureau of Astronomy reported, ‘Winter offi- 
cial Han Xianfu¥ 4% 7 completed the casting of the cop- 
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Fig. 5. Results shown by colored circles, with radii 
corresponding to calculated error. Positions of the north pole 
at different epochs are shown in blue. 
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per celestial globe.' The decree mandated its transfer to 
the Longtu Pavilion and instructed Xianfu to mentor 
selected students in the craft." The same historical source 
further details, "In the third year of Xiangfu, on the fifth 
day of the seventh month, Winter official Han Xianfu of 
the Bureau of Astronomy reported the positions of the 
outer stars relative to the polar degrees. The results were 


o=0.250 Year: 1011 Diff:0.04 
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Fig. 6. The optimal candidate time point is 1011 CE. o = 0.25 
(gudu). 


compiled into the volume 'Chongwen Mu Xiu Yao Du 
Fen Yu Ming Lu." This information is also referenced 
under the entry 'Xiangfu Stellar Degrees' in the Yuhai Æ 
305], 

According to historical records, Han Xianfu played a 
significant role in the production of the Copper Celestial 
Globe at the Zhidao Bureau in the first year of Zhidao 
(995). Engaged in crafting this celestial globe, he meticu- 
lously documented production methods, as recorded in the 
Song Shi - Tianwen Zhil!®l, "The Copper Seasonal Instru- 
ment was crafted by Han Xianfu based on the methods 
passed down from the Tang Dynasty astronomer Li Chun- 
feng4= 1% JX and the monk Yixing f$ — J. Subsequently, 
he compiled these details into a ten-volume book and pre- 
served it in the imperial library." Han Xianfu passed 
away in the sixth year of Dazhong Xiangfu, and his stu- 
dent, the Song Dynasty astronomer Yang Weide, was 
entrusted with compiling the Jingyou Qianxiang Xinshu 
Hz R Hr P (New Book of Jingyou and Qianxiang). Dur- 
ing the writing of this book, Yang Weide and his col- 
leagues conducted a series of celestial observations. 

It can be inferred that in approximately 1010 CE, the 
Winter official Han Xianfu supervised the creation of the 
Xiangfu Longtu Pavilion Copper Celestial Globe, follow- 
ing Emperor Zhenzong's instruction to teach and leave 
behind a volume of observational records. His student, 
Yang Weide, carried on Han Xianfu's legacy and orga- 
nized additional stellar observations. The data recorded in 
extant historical documents likely represent records of 
these observations, partially transmitted through subse- 
quent observations over several decades. 


4.2. Yuan Dynasty star catalog 


The National Library of China houses a Ming 
Dynasty (1368-1644 CE) copy of Tianwen Huichao K 
X Æ $b (Collection of Astronomical Material), compris- 
ing eleven sections in six volumesl!7], The third volume, 
titled Sanyuan lieshe ruxiu quji ji = 1E 7) & A Th KARE 
(Collection of Coordinates of Stars in the Three Yuan and 
Twenty-Eight Xiu), contains detailed observational records 
of approximately 740 stars. 

After its discovery and organization in the 1980s, sev- 
eral scholars undertook research on this text. In 1986, 
Chen Ying analyzed the data and, as a result of Guo Shou- 
jing’s #8 SF k stellar observations, concluded that the star 
catalog’s observations were made in the Yuan Dynasty 
(1279-1368 CE) in 1280 CENSI. In the 1990s, Sun 
Xiaochun conducted meticulous data analysis, determin- 
ing the actual observation year to be approximately 1380 
CE, during the early Ming Dynasty’s Hongwu period!]. 
He published the identification results of 678 stars. Pan 
Nai also studied and analyzed the data, reaching a conclu- 
sion similar to Chen’s, dating it to approximately 1280 
CEU, 

In this section, calculations use the data of the 669 
stars published by Sun Xiaochun. Error ranges for declina- 
tion are sampled using a Gaussian distribution. Table 2 


presents the calculated years and their different error 
radii, along with the corresponding distances (Paisti) 
between the calculated and actual north pole positions for 
each year. Fig. 7 provides visual representations at vari- 
ous error radii, while Fig. 8 shows computational results 
for the optimal case with o = 0.1 (gudu). Consequently, it 
can be reasonably inferred that the epoch is approxi- 
mately 1355 CE. 


Table 2. Yuan Dynasty star catalog (Tianwen Huicaho catalog) 
processing results 


QJD error (gudu) Epoch Paist (deg) 
0.050 1353 0.20 
0.100 1355 0.20 
0.125 1357 0.20 
0.200 1361 0.20 
0.250 1363 0.20 
0.333 1366 0.20 
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Fig. 7. Results, calculated with different error radii, are 
represented by colored circles. The positions of the north pole 
at different epochs are shown in blue. 
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The year 1355 CE corresponds to the fifteenth year 
of the Yuan Dynasty’s Zhizheng # JE era, marking the lat- 
ter part of the Yuan Dynasty. Judging from the data 
records, this set is exceptionally rare. The data is pre- 
sented in the form of star charts, with the positional data 
of stars meticulously recorded alongside them. Further- 
more, the data are organized into a percentage-based sys- 
tem, which was uncommon for that era. It can be specu- 
lated that this represents a meticulous stellar observation 
conducted by officials of the Bureau of Astronomy 
toward the end of the Yuan Dynasty, utilizing the celes- 
tial globe created by Guo Shoujing 3 “f°. These observa- 
tions were preserved in chart form and subsequently 
copied during the Ming Dynasty. 


5. CONCLUSIONS 


Based on the calculations and analysis of the two sets 
of data above, an initial evaluation of the Generalized 
Hough Transform algorithm for determining epochs can 
be provided. This method proves effective in estimating 
the observation epochs of ancient Chinese astronomical 
records and demonstrates greater accuracy compared with 
previous research efforts. The advantages of using the 
Hough transform are its independence from right ascen- 
sion data and its ability to handle erroneous data without 
specific adjustments. Through extensive error sampling, 
the impact of incorrect data can be mitigated. Simultane- 
ously, in the calculation process, various parameters 
related to stellar proper motion are taken into account, fur- 
ther enhancing the accuracy of epoch determination. 

In the calculation process of this study, observational 
data from the Song Dynasty show better data clustering 
under a condition of o = 0.25 gudu, while in the Yuan 
Dynasty, under a condition of o = 0.1 gudu, the data clus- 
tering surpasses that of the Song Dynasty. This suggests 
that from the Song Dynasty to the Yuan Dynasty, with 
the advancement of instruments, the precision of the data 
has significantly improved. 

The method employed in this study has yielded favor- 
able results on the Song-Yuan Star Catalog and offers a 
new perspective for analyzing ancient astronomical cata- 
log data in the future. 
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